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1. Introduction 

Nowadays it is well established both theoretical and experimentally that octupole 
deformation is responsible for some unusual properties of the low energy spectrum 
of several isotopes of radium and thorium as are the presence of low lying 1~ states 
in even-even nuclei pQ. As those nuclei are usually quadrupole deformed, the next 
member of the negative parity rotational band, namely the 3~ is also observed. This 
3~ is connected by strong B(E3, 3~ — > + ) transition probabilities to the ground 
state. Another typical feature of octupole deformation is the appearance of alternating 
parity rotational bands. Although there are several examples of parity alternating 
rotational bands at low spins [TJ [2] the parity alternation usually appears at high 
spin as a consequence of the stabilizing effect of angular momentum on the octupole 
character of the system [3j 0] (moments of inertia increase with the octupole moment 
and therefore configurations with higher octupole moments are favored as angular 
momentum increases). 

The appearance of octupole effects is strongly linked to the underlying single 
particle spectrum [5]. The reason is that octupolarity is enhanced when, in a given 
major shell, the intruder interacts (via a particle-hole excitation) with a standard 
parity orbital with three units less of angular momentum (3h is the amount of angular 
momentum carried out by the octupole operator). This happens between the J15/2 
and <?g/2, the 213/2 and f-j/2 an d the hn/2 and d 5 / 2 orbitals. For those regions where 
both protons and neutrons feel a strong octupole interaction is where octupole related 
effects are expected. For example, in the region around |g 4 Ra the proton's Fermi 
level is close to the f-j^ orbital that interacts strongly with the empty 113/2 nearby. 
Besides, the neutron's Fermi level is near the gg/2 orbital that strongly interacts 
through the octupole interaction with the J15/2 nearby. This combined tendency 
of both protons and neutrons towards octupole deformation is responsible for the 
observed features of negative parity states in the region around 224 Ra (see pQ for 
relevant bibliography). An interesting question is whether the octupole deformation 
effects persists when we increase proton number and move away from the ideal case 
of Ra. For neutrons we already know the answer as octupole effects are only strong 
for the radium isotopes with mass numbers in the limited range 222-228 (N=134- 
140). However, neutrons are in a different major shell in the actinides and the 
dependency with proton number is far from being the same. To answer this question 
we have carried out a series of calculations in several isotopes of the even-even 
actinides U, Pu, Cm, Cf, Fm and No exploring the possibility of octupole deformed 
ground states. To this end, we have computed potential energy curves (PECs) 
as a function of the axially symmetric octupole moment (K=0 bands) within the 
constrained Hartree-Fock-Bogoliubov (HFB) approximation [6j [7]. The PECs together 
with the collective octupole inertia are used as input for a one-dimensional collective 
Hamiltonian calculation that allows the evaluation of the energy of the lowest lying 1~ 
and 3~ states as well as B(E1) and B(E3) transition probabilities to the + ground 
state [SJ Several parametrizations of the Gogny [TU1 E] and BCP [T2l H31 IT3I [T3] 
energy density functionals (EDF) are used. The reason for the consideration of this 
variety of interactions is its empirical and phenomenological nature. The interactions 
are fitted to bulk properties of nuclei (binding energies, radii, etc) and not to the 
spectroscopic quantities considered in this paper. Therefore, a variation of their values 
with the interaction considered is to be expected. The range of variation with different 
interactions can be taken as an estimator of the error attributable to the uncertainties 
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in the nuclear interaction. 

All the nuclei considered in this work have a deep prolate minimum along the 
quadrupole degree of freedom. As a consequence the coupling between the quadrupole 
and octupole degrees of freedom is expected to be weak and therefore the octupole 
degree of freedom alone is going to play a role in the properties of negative parity states 
[T6] . There are other regions of the periodic table |17[ [T8] [19] where the quadrupole- 
octupole coupling as well as triaxiallity effects are relevant. 

Recently a theoretical survey of K=0 octupole states for even-even nuclei has 
been published [10] ■ In this survey, the Generator Coordinate Method (GCM) with 
the octupole degree of freedom has been used to compute octupole properties of over 
800 even-even nuclei including some of the nuclei discussed in this paper. However, a 
large fraction of the nuclei studied here are quite close to proton's drip line and were 
not considered in [20J. 

2. Theoretical tools 

The description of octupole properties will be carried out in the mean field framework 
[6] with the realistic Gogny [TOJ [TT] and BCP [El ES El E5] EDFs. As pairing 
correlations play a relevant role in the low energy nuclear dynamics the full HFB 
approximation will be considered. In this scheme, quasiparticle annihilation and 
creation operators are defined as linear combinations of a conveniently chosen single 
particle basis. In our case, we have used a Harmonic Oscillator (HO) basis which 
has been chosen big enough as to warrant convergence of the physical quantities 
(energies, transition probabilities, etc) with the basis size. To solve the constrained 
HFB equations an approximate second order gradient method [ST], based on the 
parametrization of the energy in terms of the Thouless expansion of the most general 
HFB wave functions and using two quasiparticle energies as precondititioner, has been 
used. Axial symmetry has been preserved in the calculation suggesting the use of an 
axially symmetric HO basis made up of the tensor product of two dimensional HO wave 
functions times one-dimensional HO ones. Along with the octupole moment constraint 
associated to the multipole operator Q% = r^Y^Q and used to generate the potential 
energy curves (PEC's), we have included a constraint on the center of mass of the 
nucleus (i.e. the mean value of r 1 Yio has been set to zero) to prevent the coupling to 
the spurious center of mass motion. As a consequence of the axial symmetry imposed 
in the HFB wave functions the remaining components fi ^ of the corresponding 
multipole operators r 3 Y^^ and rY\^ have a zero mean value by construction. 

The information given by mean field theories is restricted to the energy and shape 
of the -generally- deformed ground state. To restore the parity symmetry broken by the 
solutions of the mean field approximation and in order to describe the dynamics of the 
collective excited states of octupole character it is mandatory to go beyond the mean 
field approximation. With this in mind, the octupole degree of freedom Q3 = (i/j\Q3\ip) 
(where is the HFB intrinsic wave function) has been used to build up a collective 
Hamiltonian based on the GCM and the Gaussian Overlap Approximation (GOA) 
[221 1231 [24] . In this method, the GOA is used to reduce the Hill- Wheeler equation 
of the GCM to a Schrodinger equation for the collective wave function, the so-called 
Collective Schrodinger Equation (CSE) 

Ucoll<j>a{Qz) — £a4>a(Q?,), (1) 
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+ V(Q 3 )-e (Q 3 ). (3) 



In this expression G{Q 3 ) is the metric, B(Q 3 ) is the mass parameter associated 
with the collective motion along Q 3 , V(Q 3 ) is the collective potential given by the 
HFB energy V(Q 3 ) = (xp(Q 3 )\H\ip(Q 3 )) and e (Q 3 ) is the Zero Point Energy (ZPE) 
correction. The eigenfunctions <j) a {Q 3 ) of equation ([1]) have to be normalized to one 
with the metric G(Q 3 ) 

J dQ 3 y/G(Qd<l>* a (Q3)MQs) = <W (4) 

in order to preserve the hermiticity of ticoii ■ 

It has to be mentioned that a collective Schrodinger equation can also be obtained 
from the Adiabatic Time Dependent Hartree-Fock (ATDHF) theory [25l EH [27] after 
quantization of the semi-classical Hamiltonian for the slow moving collective degrees of 
freedom. The collective Hamiltonian obtained in this way has the same functional form 
as the GCM+GOA one, but the expression of the collective parameters is different. 
Later on we will discuss how to choose these collective parameters. 

An interesting property of % co u is its invariance under the exchange Q 3 — > — Q 3 
that allows the classification of its eigenfunctions, (f> a (Q 3 ), according to their parity 
under the Q 3 —> —Q 3 exchange. It is easy to see that the parity of the collective wave 
function under the Q 3 —> —Q 3 exchange corresponds to the spatial parity operation in 
the correlated wave function built up from cj> a . The inclusion of octupole correlations 
immediately restores the parity symmetry lost at the mean field level [8l[28|. Therefore, 
the solution of the CSE equation (fTJ allows the calculation of the + — 1~(3~) energy 
splitting as well as the B(E1) and B(E3) transition probabilities connecting them. 
At this point it has to be pointed out that in the present framework where only time 
reversal invariant wave functions are considered it is only possible to describe excited 
states with average angular momentum zero. To deal with genuine or 3 _ states 
a projection onto good angular momentum should be performed, which is out of the 
scope of the present work. Here we will assume that the cranking rotational energy 
of the Instate is much smaller than the excitation energy of the negative parity band 
head and therefore can be safely disregarded. With this approximation in mind, the 
reduced transition probabilities from the lowest 1~ and 3~ states to the + ground 
state can be computed within the Rotational Model approximation as 

B(EX,I f -> h) = e 2 (I i KX0\I f K) 2 \(cp i \r x Y Xfi \ t p f }\ 2 , (5) 

where \ipi) and \iff) are correlated wave functions obtained in the spirit of the GCM 
from the collective wave functions 4> a {Q 3 ). The above formula can be reduced to an 
expression involving those collective wave functions (f) a (Q 3 ) by means of the GOA [30J. 
The final result for K = bands reads 

|A)|0o+>coll| 2 (6) 



|Q 3 o(PROT)|0 o+ ) CO ll| 2 , (7) 



B(£1,1-^O+) = ^K0 O - 
for the El electric transition and 

B(E3, 3-^0+) = ^-\(<t> - 
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for the E3 one. In the above formulas we have introduced the collective matrix element 
of an operator O as 

(0 o -|6|0 o+ )coll = j dQ 3 G 1 ' 2 <&-(Q a )0(Q a )fa(Q 3 ) 

where 0(Q 3 ) = (^(Q a )\0\ip{Q a )) . In the formula © Dq is the dipole moment 
operator whose mean value is defined as the difference between the center of mass 
of protons and neutrons 

N Z 

D (Q 3 ) = -jmQ 3 )\z P roMQ3)) - ^(QsjlwI^Qs)). (8) 

Finally, Q 3 o(PROT) is the part of the octupole operator acting on proton's space. In 
the recent survey of |20| it was pointed out the inadequacy of the rotational formula 
for the calculation of El and E3 transition probabilities for near spherical nuclei (see 
|29) for a beautiful example of how angular momentum projection solves a related 
problem near 208 Pb ). This is not a problem in the present calculation as the nuclei 
considered are well deformed and the rotational formula works well there. 

To carry out the collective calculations it is necessary to specify the collective 
parameters G(Q 3 ), B{Q 3 ) and eo(Q3) appearing in the definition of Hcoii equation 
©. As it was mentioned before, there are two sets of parameters coming from the 
GCM+GOA and the ATDHF derivation of the collective Hamiltonian. The set of 
parameters used in this calculation is an admixture of the two and it is known as 
the ATDHF+ZPE set. It includes the mass parameter B(Q 3 ) coming out from the 
semi-classical Hamiltonian of the ATDHF theory, the metric of the GCM+GOA and 
the ZPE correction computed with the GCM+GOA formula but using the ATDHF 
mass instead, i.e. 

eo(Q 3 ) = Ig(Q 3 )B(Q 3 ) a ^ dhf . (9) 

This set of parameters was devised to put together the advantages of the ATDHF set 
(time-odd components included in the mass term) and the ones of the GCM+GOA 
(ZPE correction). This method can be somewhat justified in the context of the 
extended Generator Coordinate method |24[ I3T] and has been extensively used |8] . 
In mean field calculations it is customary to include in the ZPE correction the 
rotational energy. We have tested that the inclusion of the rotational energy correction 
(as computed in to take into account effects near sphericity and the approximate 
calculation of the Yoccoz moment of inertia) does not modify substantially the results, 
as could be anticipated due to the weak coupling of the octupole and quadrupole 
degrees of freedom. 

The calculation of the collective parameters involves the inversion of the HFB 
stability matrix which is closely related to the matrix of the RPA equation. At present, 
this is a formidable task and approximations are needed. The approximation used in 
this paper - called "cranking approximation" |32[l33] - neglects the off-diagonal terms 
of the stability matrix allowing to invert it analytically but at the cost of including 
the two body interaction only through the mean field. Although this approximation 
has been extensively used in the literature for the calculation of collective masses and 
moments of inertia (see for instance |34[ 1111 135| ) its validity has not been properly 
established. Using the "cranking approximation", the ATDHF+ZPE parameters are 
given by 

, . _ M_ 2 (Q 3 ) _ M_ 3 (Q 3 ) . 

G(Qs) - 2M^(g 3 )' B(Q3) - wm' m 
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where the quantities M_ n (Q 3 ) (n = 1,2, 3 ) are defined as 



M_ n (Q 3 ) = 



E 

k,l 



\(Q30)kl 



(E k +E[) 



(11) 



In the above expression, E). are the quasi-particle energies and (Q3o)t.; are the matrix 
elements of the 20 part |6j of the octupole operator Q30 in the quasi-particle basis of 
the HFB wave function \4){Q^j). This form of the collective mass is usually referred 
to in the literature as the Belyaev-Inglis inertia [5J. 



2.1. Interactions 

In this paper we have used the three modern parametrizations of the Gogny functional 
[TU] which are available in the market, namely Gogny D1S [II], DIN [37] and DIM [55] . 
There are mainly two reasons to proceed in this way. First, to gain confidence on the 
independence of our results concerning tiny details of the interaction/parametrization 
used. In this respect, we have to comment that although the three Gogny's 
parametrizations used belong to the same functional form of the interaction, their 
different behavior regarding nuclear matter, pairing properties, binding energies or 
radii suggest their non-equivalence. The other reason to do the calculations with the 
DIN and DIM parametrizations is to check their performance in describing octupole 
deformation phenomena as a partial test to asses their ability to describe low energy 
nuclear phenomena (see for instance for recent studies with DIM). In addition 

to the Gogny force, we have used the recently proposed BCP [12] energy density 
functional with the BCP1 parameter set [TJ]. This functional is based on a local 
density approximation to a realistic nuclear matter equation of state supplemented 
by a finite range surface term. The spin-orbit and Coulomb interactions have the 
standard form as used in the Gogny force. The pairing force used is a zero range, 
density dependent force that was fitted to reproduce the pairing gap in nuclear matter 
given by the Gogny functional. The parameters have been fitted to reproduce binding 
energies and radii of selected spherical nuclei and the results show a reasonable level 
of agreement comparable or even better than the results provided by D1S. Several 
spectroscopic quantities have been computed in several regions of the periodic table 
and the results are in close agreement with the experimental values and the theoretical 
results of D1S [13] H5] . Octupole properties of Ra, Th and Ba isotopes have been 
analyzed in the past with this functional [14] and the good agreement obtained at a 
qualitative level in the comparison with the experiment and Gogny force calculations 
is a clear indication of the suitability of this EDF for the exploration of the octupole 
degree of freedom in the heavy actinides. 



3. Results 



In the eight isotopic chains considered we have followed the same procedure, first 
the octupole constrained HFB wave functions have been generated along the lines 
described in the previous section. Next the collective mass and zero point energies are 
computed with the mean field wave functions and together with the PEC itself they are 
used to build the Collective Schrodinger Hamiltonian. From the lowest lying solutions 
of the CSE we obtain the excitation energy of the first negative parity excited state 
(a \- for these deformed nuclei) and the B(E1, 1~ — > 0+) and B(E3,3~ — > 0+) 
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Figure 1. On the left hand side panels, the PEC's as a function of the octupole 
moment Q3 (in thousands of Fermi cube or 6 3 ' 2 ) for the even-even isotopes of 
uranium with A in the range between A = 220 and A = 240. Results for the D1S 
(full lines), DIN (dashed lines), DIM (dashed dotted lines) parametrizations of 
the Gogny and the BCP1 (dotted lines) functionals are presented. On the right 
hand side panels the particle-particle correlations energies for both protons (thick 
curves) and neutrons (thin curves) are given as a function of the octupole moment 
in b 3 /' 2 = 10 3 fm 3 . 

transition probabilities (the 3~ state in the B(E3) transition probability refers to 
the second member of the rotational band built on top of the l~band head). The 
whole procedure will be thoroughly described for the case of the uranium isotope. 
For the other isotopic chains only the most relevant results, namely the PES and the 
observables obtained after solving the CSE will be presented. 

3.1. Low excitation energy properties of the uranium isotopes 

In fiugre [T] we display, as a function of the octupole moment Q3 and the four 
interactions used, the PECs (left hand side panels) and the particle-particle correlation 
energies (right hand side panels). The particle-particle correlation energy is defined 
as E pp = iTrA/i, and it is given in terms of the pairing field A and the pairing 
tensor n. This quantity gives a rough idea of the amount of pairing correlations in 
the system. It can also be used as an indicator of where the single particle level 
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density around the Fermi surface is high or low as high level densities are associated 
to strong pairing correlations. This quantity is also correlated with the pairing gap 
that represents the energy of the lowest two quasiparticle excitation and therefore it 
is closely related to the collective inertias to be discussed below. Coming back to the 
figure, the PECs obtained with the four interactions considered are very similar. For 
the particle-particle correlation energies the same similarity is obtained for the three 
parametrization of Gogny, but for the BCP functional the particle-particle correlation 
energies are significantly smaller. Of the eleven isotopes considered, five ( 222 - 230 U) 
show an octupole deformed minimum with a depth never exceeding the 1.5 MeV. The 
octupole moment of the minima changes with the neutron number in a non-trivial 
way. The remaining nuclei with a minimum at Q3 = show a parabolic behavior 
centered around the minimum with a curvature that is more or less independent 
of the interaction considered. There seems to be a certain correlation between the 
appearance of octupole deformed minima and the minima of the neutron pairing 
energies of the right hand side panels. This is something to be expected as the presence 
of deformed minima is usually related (Jahn- Teller effect) to regions of low level density 
of the underlying single particle orbitals and these low level density regions usually 
correspond to small pairing correlations. 

As the octupole moment increases, the nuclear shape changes and the values of 
other multipole moments including the quadrupole, hexadecapole and dipole moments 
change. In the left hand side panels of figure [5] the fa and /?4 deformation parameters 
are represented as a function of Q3. For most of the U isotopes considered the 
variation with Q3 is small and it is only in the lightest nuclei 220 U and 222 U where 
the changes are noticeable. The nucleus 220 U is spherical at Q3 = and develops a 
sizable quadrupole deformation with increasing octupole moment. On the left hand 
side panels of this figure the dipole moment of equation © is plotted. The dipole 
moment is zero by construction at Q3 = and its behavior with the octupole moment 
is determined by the occupancies of specific single particle orbitals around the Fermi 
level. High K orbitals give the largest contributions to Dq and therefore, depending on 
their occupancy for protons or neutrons as a function of Q3, we can observe different 
behaviors. For the lightest isotopes the dipole moment increases almost linearly with 
Q3. The slope decreases with neutron number to the point that for 236 U the value of 
Do is almost zero for a wide range of octupole moments as a consequence of competing 
neutron and protons contributions. Given the intimate connection between the dipole 
moment and the B{E\) transition probabilities we conclude that the situation found 
in 236 U will lead to a quenched B(E\) value as compared to the ones of other isotopes 
of the same nucleus. 

In figure |3] the collective inertia B(Qs) of equation (fT0|) is depicted as a function 
of Q3. It plays a central role in the outcome of the collective Schrodingcr Hamiltonian 
calculations mentioned in the previous section. The collective inertia is roughly 
speaking inversely proportional to the pairing correlations (the pairing gap to be 
more specific) and therefore the similarity between the collective masses obtained 
for D1S, DIN and DIM indicate that a subtle cancellation between proton's and 
neutron's contributions has to take place (see the particle-particle correlation energies 
in figure [lj. It is also worth to notice that the peaks observed in the B(Qs) plots are 
related to regions of low pairing correlations. On the other hand, the BCP1 inertia 
is systematically higher than the ones of the different parametrizations of the Gogny 
functional. The difference amounts to factors of two or even more in some specific 
cases. This will be responsible for the differences observed in the solution of the CSE. 
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Figure 2. On the left hand side panels, the 02 (thick lines) and 04 (thin lines) 
deformation parameters are plotted as a function of the octupole moment Q3 (in 
£,3/2 = 10 3 fm 3 ) for the even-even isotopes of uranium with A in the range between 
A = 220 and A = 240. Results are given for the D1S (full lines), DIN (dashed 
lines), DIM (dashed dotted lines) and BCP1 (dotted lines) functionals. On the 
right hand side panels the dipole moment is shown as a function of Q3. 

On the right hand side panels of figure [3] the zero point energy correction e(Q3) of 
equation is depicted for the isotopes of uranium. The e(Qz) values are correlated 
with the inverse of the collective inertia B(Q^) as can easily be noticed in the figure 
and its behavior with Q3 can modify somehow the topology of the PES but not at the 
level of changing significantly the physical results. On the average the e(Qz) values 
are of the order of 1.5 MeV which is consistent with the correlation energies for parity 
symmetry restoration and octupole dynamics obtained in [20J. This energy together 
with the energy gain obtained by the breaking of reflection symmetry at the mean field 
level (of the order of one MeV at most) should be added to the nucleus binding energy 
but the relatively constant values as a function of mass number indicate that the 
impact of such addition in separation energies is not expected to be significant. After 
solving the one-dimensional collective Schrodinger equation we obtain the excitation 
energies of the 1~ states and the B(E1, 1~ — > + ) and B(E3,3~ — > + ) transition 
probabilities. The results obtained with Gogny D1S, DIN and DIM as well as for the 
BCP1 energy density functional for the U isotopes are depicted along with available 
experimental values [JTjin the left hand side panels of figure 01 For all the results 
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Figure 3. On the left hand side panels the octupole collective inertia parameter 
B{Q'i), which is one of the relevant ingredients in the one-dimensional collective 
Schrodinger Hamiltonian (see text for details), is shown as a function of Qg for 
the considered U isotopes. On the right hand side panels the zero point energy 
(ZPE) correction of equation J9} is represented as a function of Q3. 



considered we observe that the three parametrizations of the Gogny EDF produce the 
same isotopic trend being the differences not very relevant from a physical point of 
view. Obviously, the differences are larger for the most sensitive quantities, namely 
the transition probabilities. The 1~ excitation energies are very small for the 222 U to 
230 U isotopes. This is the expected behavior as those nuclei show permanent octupole 
deformation at the HFB level (see figure[T]). For the other U isotopes, the HFB energy 
shows a minimum at Q3 = with a parabolic behavior characteristic of octupole 
vibrational states which lie higher in energy, as observed in the energies of the negative 
parity excited states of figure @] The comparison with experimental data is reasonable 
as we are able to reproduce the increasing octupolarity as neutron number decreases. 
Concerning the B(E1, 1~ — > + ) transition probabilities we observe a pronounced 
minimum around A — 236 that is a direct consequence of the behavior of the dipole 
moment with Q3 observed in figure [5J Unfortunately, there is no experimental data 
available for these isotopes concerning the B(E1) transition probabilities and the 
associated dipole moments [32]. The B(E3,3~ — > + ) transition probabilities show a 
marked maximum around A — 226 that is correlated to the minimum in the energies 
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of the 1~ states. For the heavier isotopes the computed B(E3) compare well with 
experimental data. On the other hand, the BCP1 results compare reasonably well 
with the Gogny ones except for the 228 - 232 U isotopes where the octupole correlations 
are weaker and therefore the 1 _ excitation energies are higher and the B{EZ) transition 
probabilities lower. Also the minimum observed in the B(E1) as a function of neutron 
number is shifted two units of neutron number when compared with the Gogny results. 
In the right hand side panels, a comparison is made between the theoretical results 
obtained in the present framework and the ones of the GCM with the octupole 
degree of freedom as generating coordinate (GCM-Q3) obtained for Gogny D1S in 
|20) . In both calculations the underlying mean field is the same and therefore it 
should not be surprising the agreement between the results of both calculations. The 
differences observed in the excitation energies and B(E3) transition probabilities can 
be attributed to the different collective masses used in each of the calculations 0. 
However, the magnitude of the differences is compatible with the magnitude of the 
differences attributed to the use of different parametrizations and/or functionals in the 
calculations. The message that the comparison of all these theoretical predictions is 
telling us is that the behavior of the physical quantities as a function of mass number 
is independent of the kind of calculation and interaction and therefore can be taken as 
a strong prediction. On the other hand, the values obtained for the physical quantities 
oscillate typically in a range of plus or minus 40 % and therefore their values should 
be taken just as an indication of the order of magnitude to be expected. 

3.2. Plutonium isotopes 

Results for the Plutonium (Z=94) isotopes in the mass range 222-242 (corresponding to 
neutron numbers 128 - 148 ) are discussed in this section. For this and the remaining 
isotopic chains to be discussed, we will focus on physical quantities like excitation 
energies and transition probabilities and the only "intrinsic" quantity to be shown is 
the PES as a function of the octupole moment to quantify the gain in binding energy 
associated to octupole correlations. In figure [5] we have plotted on the left hand 
side the PES as a function of the octupole moment for the Pu isotopes considered. 
Well deformed octupole minima are obtained for the isotopes 224 Pu, 226 Pu and 228 Pu 
with mean field octupole correlation energies as large as 1 MeV in 226 Pu. The nuclei 
230-234p u snow ver y shallow minima which are located at Q30 = in the heaviest 
of the three. For the other isotopes, the minima are reflection symmetric (Q30 = 0) 
showing a parabolic behavior as a function of Q30 with a rather large curvature. 
On the right hand side the 1~ excitation energy as well as the B{E1) and B(E3) 
transition probabilities (in Weisskopf units) are plotted as a function of mass number. 
We observe that a range of isotopes starting with 22 Pu and ending with 232 Pu show 
octupole deformation in their ground state. The mean field octupole correlation energy 
is of the order of 1 MeV in the lighter isotopes and gets reduced to a few hundred 
keV for 232 Pu. As a consequence, the 1~ excitation energy for those five isotopes is of 
the order of 0.5 MeV and much smaller than in the neighboring isotopes. As a result 
of the enhanced octupolarity in these isotopes the B(E3) transition probabilities are 
stronger reaching values of up to 60 W.u. The B{E1) transition probabilities show a 
rather constant behavior for the octupole deformed isotopes and show a dip for 238 Pu 

\ In the GCM-Q3 calculation there is no explicit collective mass as the full Hill- Wheeler equation 
is solved but, if a local approximation to the non-local HW equation is considered the concept of 
collective mass can be recovered, see Ref [6] for a discussion. 




Figure 4. Left hand side panels: Results of the ID collective Schrodinger 
Hamiltonian calculations for the three parametrizations of Gogny (D1S, DIN and 
DIM) as well as for the BCP1 energy density functional. In the lowest panel the 
energies of the 1 _ states as a function of the mass number for the uranium isotopes 
are depicted. Theoretical results (lines) are plotted along with experimental data 
(dots). In the middle panel the B(E1, 1~ — > + ) transition probabilities in W.u. 
are given. Finally, in the upper panel the B(E3, 3~ — ¥ 0) are also given, for the 
different isotopes and different theoretical prediction, in W.u. In the right hand 
side panels a comparison is made between the results obtained in the present 
approach and the ones obtained with Gogny D1S in the framework of the GCM 
for the octupole degree of freedom. 



in the case of D1S and 240 Pu in the case of both DIN and DIM. 
3.3. Curium results 

Results for the Curium (Z=96) isotopes with neutron numbers in between N=126 and 
N=146 are discussed in this section. In figure [5] we have plotted on the left hand 
side the PES as a function of the octupole moment. Octupole deformed minima are 
obtained for the nuclei 22e ~ 230 (j m with energy gains of the order of a few hundreds of 
keV being 226 Cm the nucleus with the deepest octupole deformed well. For 232_234 Cm 
the minima are at Q30 = but the minima are in the two cases rather shallow. For the 
other nuclei the minima is at Q30 = and the PES shows a parabolic behavior with a 
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Figure 5. On the left hand side panels, the PES, as a function of Q3, are shown 
for the considered Pu isotopes. In the right hand side panels, the 1~ excitation 
energies (in MeV) as well as the B(E1) and B(E3) transition probabilities (in 
W.u.) are plotted as a function of the mass number A. Theoretical results are 
depicted as lines of various types (see legends) and the experimental data, taken 
from 1411 . are represented by bullets. 



rather large curvature. On the right hand side panels, the 1~ excitation energy as well 
as the B(E1) and B(E3) transition probabilities (in Weisskopf units) obtained after 
solving the ID CSE are plotted as a function of mass number. The results resemble 
at a qualitative level the ones for the Pu isotopes indicating the relatively low impact 
on octupole correlations of the addition of two protons. Strong octupole correlations, 
with low lying l - states are anticipated to show up in the mass range 226-232. In that 
range the B(E3) transition probabilities are rather collective and are around 50 W.u. 
For the nucleus 234 Cm the results show a dependency with the interaction, pointing 
to a characterization of this nucleus as transitional in terms of octupole correlations. 



3.4- Californium 

Results for the Californium (Z=98) isotopes with neutron numbers in the 124 -144 
range are presented and discussed in this section. In figure [7] we have plotted on the 
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left hand side the PES as a function of the octupole moment. On the right hand side 
the 1~ excitation energy as well as the B(E1) and B(E3) transition probabilities (in 
Weisskopf units) are plotted as a function of mass number. The qualitative and even 
quantitative similarity with the Cm results presented in the previous section is striking 
and indicates that the addition of the two protons has little influence in the octupole 
deformation properties. Strong octupole correlation effects are predicted in the mass 
number range 228 - 236 and a severe quenching of B(E\) transition probability is 
anticipated for the nuclei around 224 Cf. At this point it has to be mentioned that 
proton's drip line is predicted to be at 228 Cf in calculations with Gogny DIS [43J. 
Therefore it is rather unlikely that the octupole effects expected in the range 228-236 
will ever be accessible experimentally. 



3.5. Fermium and Nobelium 

For the Fermium (Z=100) and Nobelium (Z=102) isotopes calculations with neutron 
numbers in the range N=122 and N=142 have been carried out. Although octupole 
deformation is present in our calculations for the ground state of 23 °- 232 Fm (N = 
130 - 132) and 2 32-234 No ^ N _ 13Q _ 132 ) we have decided not to present the results 
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Figure 7. Same as figure [5] but for the californium isotopes 



here as those nuclei are right at the proton's drip line, according to calculations with 
the Gogny D1S force (see [43] for details), and therefore it is rather unlikely that they 
will ever be created in the laboratory. 

4. Conclusions 

We have performed systematic mean field calculations for a relevant set of even-even 
actinides to explore the impact of octupole correlations in the low energy nuclear 
spectrum. We observe strong octupole correlations in some uranium, plutonium, 
curium and californium isotopes with a strength that decreases with proton number. 
The neutron number of the isotopes with strong octupole correlations become closer to 
the proton's drip line as proton number increases. As a consequence, the fermium and 
nobelium isotopes showing up noticeable octupole correlations are beyond the drip 
line. Excitation energies for the l~state as low as 300 keV and B(E3) transition 
probabilities of 60 to 70 W.u. are predicted in some uranium isotopes, making 
them good candidates for experimental studies aiming to extend the region of known 
nuclei where octupole correlation effects are important. On the theoretical side, the 
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similarity among the results obtained with three different parametrizations of the 
Gogny functional give credit to them as suited for an equivalent description of octupole 
correlation effects. The similitudes between the Gogny and BCP results indicate 
that the octupole correlation effects obtained are genuine and not an artifact of the 
phenomenological interactions used. 
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